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An asymptotic description of the laminar flow pattern generated by a point source of 
heat is derived. Existing solutions of the boundary layer equations are shown to be valid 
in a limited domain where buoyancy forces are dominant. These solutions are supplemented 
by new results for the buoyant plume and are matched asymptotically to solutions valid 
everywhere outside the plume. Composite analytical formulae for the pressure, enthalpy 
and velocity fields are obtained and applied to the computation of the streamline pattern, 
pressure, and enthalpy distribution for Prandtl numbers 0.7 and 1.0. 
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1. Introduction 

The flow induced by an isolated heat source is a prototype for many buoyancy driven 
flows. The motion is initiated because the heated gas adjacent to the source becomes less dense 
and rises. The rising hot gas column draws in cool stationary fluid from the surroundings, 
reducing the excess temperature and hence the relative buoyancy in the moving fluid. The 
resulting steady flow pattern is a long narrow vertical plume of heated fluid rising through a 
nearly quiescent medium. The intuitive picture sketched above has led most previous investi- 
gators to study the problem by applying the boundary layer approximation to the equations 
of motion appropriate to a weakly buoyant flow (the Boussinesq equations). The procedure 
results in a system of ordinary differential equations which has been solved many times. Thus, 
Schuh [1] l obtained numerical results for a fluid of Prandtl number (Pr) 0.7. i r ih [2] gave closed 
form solutions for Pr=l and Pr=2. Fujii [3] presented numerical results for Pr— 0.01, 0.7, 
and 10, in addition to the closed form solutions for Pr=l and 2. Most of Fujii's results were 
repeated independently by Brand and Lahey [4]. A useful survey of natural convection phe- 
nomena containing these and related references is given by Gebhart [5] 

The above analyses cannot describe the entire flow pattern, however. The fluid entrained 
in the plume must come from somewhere; hence the surrounding fluid cannot be at rest. 
Moreover, the pressure cannot be constant, not even within the plume. In fact, inspection of 
the solutions to the boundary layer equations raises a question as to the sense in which boundary 
layer theory is applicable to the present problem. The difficulty arises from the fact that the 
solutions decay algebraically to their ambient values outside the plume. This behavior is very 
different from the exponential decay which obtains in most boundary layer problems, including 
the two dimensional (line source) thermal plume [3]. 

The situation is analogous to that which occurs in the study of the flow induced by a 
momentum source. The plane momentum jet (line source) was studied using the boundary 
layer equations by Schlichting [6] and Bickley [7]. An English language account of Schlich ting's 
work may be found in Schlichting [8]. The velocity profile in the plane jet decays exponentially 
as the ambient fluid is approached. The corresponding boundary layer solution for a point 
momentum source [6] is again algebraic, as in the problem under consideration. Fortunately, 
the exact solution of the Navier Stokes equations for a point source of momentum has been 
given by Landau [9]. The boundary layer solution describes the velocities (but not the pressure 
distribution) near the axis of symmetry provided the cylindrical coordinates that appear in 
the solution are reinterpreted as small angle approximations to spherical polar coordinates. 
However, the flow away from the axis is nontrivial and cannot be obtained from the boundary 
layer equations. 



1 Figures in brackets indicate the literature references at the end of this paper. 
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In the present paper an asymptotic solution for the flow induced by the point heat source 
is derived. The solution is valid sufficiently far from the source in all directions, unlike the 
solutions based on the boundary layer equations. The key elements in the solution are the 
determination of the pressure distribution in the plume and the use of parabolic coordinates. 
The introduction of parabolic coordinates is suggested by the work of Landau [9] and Goldstein 
[10] on the point momentum source and flat plate boundary layer, respectively. It permits a 
simple derivation of the boundary layer form of the equations which apply in the plume region 
above the heat source. These equations are solved analytically for any Pr of order unity. They 
lead to an expression for the pressure in the plume, a quantity not previously calculated. The 
form of the pressure distribution, together with corresponding expressions for the velocity and 
enthalpy in the plume, allow simple analytic expressions to be obtained for the velocity and 
pressure distributions in the fluid both inside and outside the plume. Finally, the enthalpy 
distribution outside the plume is determined and combined with the solution in the plume to 
produce a result valid in all directions far from the source. 

Fendell [11] treated the flow induced by a small heated sphere, and a portion of that analysis 
is closely related to the work presented here. Fendell derived the boundary layer equations for 
the velocity and enthalpy in parabolic coordinates, but assumed the pressure to be constant. 
The flow pattern outside the plume was determined using an Oseen linearization about the 
maximum plume velocity. This idea, while attractive in principle, leads to solutions which are 
qualitatively different from those obtained below. The use of the Oseen linearization in free 
convection problems appears to be more difficult than is the case in forced convection. Even 
when forced convection occurs, bouyancy can still introduce difficulties as Wesseling [12] has 
shown. 

2. The Equations of Motion 

The Boussinesq form of the Navier Stokes equations for the steady flow of a constant 
property fluid in a gravitational field may be written in vctor invariant form as: 



(H- 



(1) 



V-u=0 

-UX<o+-V(ff-ffJ + 7 m) g=-vVXa 

VXU=co. 

The first three of equations (1) are the conservation of mass, momentum, and energy, respec- 
tively, while the final equation is a definition of the vorticity w. The dependent variables are 
the fluid velocity u, the pressure p, and the enthalpy per unit mass h. The fluid properties are 
the density p, the kinematic viscosity v, and the Prandtl number Pr. The quantity g is the 
gravitational acceleration, and the subscript °° denotes the value in the undisturbed fluid far 
from the source. The boundary conditions are: 

VL=p—p m =h—h m =0; r->oo 

f Lu-nih-hJ -p- r vA-nj cIs=2tQ. (2) 

The boundary conditions (2) assume that the heat source of strength 2ttQ is located at the 
origin of coordinates r = 0. The integral is taken over an arbitrary surface which encloses the 
origin. 

The motion is symmetric about an axis parallel to the gravity vector g which passes thru 
the source. Let r denote the radial distance from the axis of symmetry and z be the distance 
above the source measured along the axis. The buoyancy dominated region lies above the 
source with a radial extent given by [3] : 

r~(z) 1/2 . 
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The motion induced everywhere else is essentially irrotational. The configuration is illustrated 
schematically in figure 1. As mentioned earlier, it is convenient to introduce parabolic coordi- 
nates defined as follows : 



£i={Vr 2 +z 2 +2} 1/2 

Let i u i 2 , be unit vectors in the direction of £ b £ 2 ; and i 3 be defined by: 

i3=i!xi 2 . 

The first of equations (1) can be used to define a stream function ^. 

u=uii+t?i 2 

1 1 d* 



(3) 



V$i 2 +fc 2 && d & (4) 

_ -1 1 d£ 



The conservation of momentum (the second of equations (1)) becomes 
du dv I 2 . , .16, . 

. . f du . dv . / . 9 , . o ,19, J 

[ o£ 2 o& p o£ 2 J 



(5) 



-M*(£- 1 )-iSK* w ) 

The terms on the left-hand side of eq (5) are the inertial and pressure contributions to the 
momentum balance while the right-hand side contains the effect of buoyancy and viscosity. 
The quantity g is the magnitude of g, and co is related to the vorticity « and the stream function 
^ as follows: 



w = l 3 co 

(Sl'+fc 2 ) 



-j 3 [ d / i a*\ a / 1 d*\\ 



Finally, the energy equation becomes: 

Q=h-h m . 



(7) 



The left-hand side of eq (7) represents the convection of energy due to the fluid motion, while 
the right-hand side accounts for the diffusion of heat. 
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£, = CONSTANT 




£2= CONSTANT 



Figure 1. Schematic drawing of coordinate geometry and plume orientation. 

The solution procedure is based on the fact that the buoyant plume drives the flow. Pre- 
vious analyses cited above show that the buoyancy dominated region lies in the domain defined 
by £i^£ 2 (see fig. 1). The appropriate dimensionless variables are: 

v _(±gQ\'\ 
V -\^K?) & 



t-\i P hy) ?1 



env) 



(8) 



pv \pv 6 h a / 
_gQ P(n) 



1/2 



6(v) 



p-p a = 



•h a ^ 



The differing factors of (4) 1/4 appearing in the definitions of f and n are chosen so the boundary 
layer equations valid for £^77 are the same as those obtained by Fujii [3]. The velocity and 
vorticity components u,v and co are obtained from equations (8) as : 



,-fM 



(M_Y 

\p V 3 hJ 



V i+ 



ldF 



Irfydr] 

II" 2 






^r 



F(n) 



\pv 3 h a J A 1 jA dri \v dv / 



"4 e) 



(9) 
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Substitution of equations (8) and (9) into eqs (5) and (7) yields, for £»*/: 

dP.nO F£/F\_ n /l\ 

Thus, letting £ approach infinity with 77 held fixed, the right-hand sides of equations (10) may 
be ignored. The first two of these equations, the i x component of the momentum equation and 
the energy equation, have been obtained many times before, but with £ and 77 replaced by their 
limiting values for z^>r: 

£~V2i 
r 

The last of equations (10), which represents the i 2 component of the momentum balance, has 
apparently never been derived previously. Its importance will become apparent below. 
The boundary conditions (2) become, to the present degree of approximation: 

P=^=G=0; „->» (11) 

P" r,dF . , 
Jo G lf= l 

The conditions at 77 = insure regularity of the solutions on the axis of symmetry. The con- 
ditions at infinity are strictly valid only if the solutions so obtained can be continued analy- 
tically, so that equations (2) can be satisfied at infinity for all values of £/t7 outside the plume. 
The integral constraint on the solution is consistent with the other boundary conditions. 
Physically, it states that the heat generated at the source is primarily convected up the buoyant 
plume. The role of thermal conduction of heat only enters the overall energy blance (the last 
of equations (2)) in terms of higher order in (£) -1 . 

3. Solution of the Boundary Layer Equations 

As noted in the introduction, there are many numerical and analytical solutions available 
to the first two of equations (10) for particular values of the Prandtl number. Rather than 
obtain additional numerical solutions, approximate analytical solutions are constructed. They 
are valid for any Prandtl number of order unity, which includes most common gases. The 
second of equations (10) may be integrated once. Applying the boundary conditions as 77 
approaches infinity and defining a new independent variable y, the first two of eqs (10) become : 



w d?F. d r d 2 Fl, r n 



HWJjf-0 02) 

v 2 =2y 
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The method of solution rests on the assumption of the following functional form for F: 

Equation (13) is suggested by the fact that the exact solutions for Pr=\ and Pr=2 have this 
form. Moreover, the boundary conditions at y=0 and the asymptotic behavior at infinity are 
satisfied for any value of Pr. Finally, eq (13) satisfies the requirement that the solutions be 
analytic in y near y= 0, a requirement imposed by the boundary conditions. Once equation (13) 
is adopted, G can be determined immediately from the second of eqs (12). 

£=<7o(l+/ty)- x 

(14) 
,_Pra_PrF(*>) 

The single additional parameter g is related to a and /3 by the overall energy balance. Using 
eqs (13) and (14) in the integral in equation (11) : 

,(i + f)"W 

The remaining parameters a and /3 are determined by requiring the first of equations (12) (the 
momentum equation) to be satisfied on the average and at y=0. This is a widely used proce- 
dure for obtaining approximate solutions in boundary layer theory L 8]. At y=0: 

d 2 F 

£(0)+2~(0)=0 

Thus: (15) 

4a/3=0 o 

The integrated i x momentum equation may be put in the form: 

fKf)>=° 

Substituting equations (13) and (14) into this expression: 
Equations (13)- (16) may be readily solved to obtain: 



a=« 



1J3+V1+24 Pr 



2 2 



1/2 



= ^yr+247y [l+VIT ^ r (17) 



p r 

ffo=-o[3+Vl+24 Pr] [l+Vl+24Pr] _1 



X=i[l+Vl+24Pr] 



Equations (13), (14) and (17) constitute the approximate solution for F and G. The solution 
is exact for Pr=l and 2. Its accuracy was tested by comparing the present results with those 
of Fujii [3] for Pr = 0.7. The comparison is shown in the following table. 
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Quantity 


Fujii [3] 


Present 


G(0) 


0.938 

7.528 
.481 


0. 950 

7.456 

.484 



The entire velocity profile was also plotted against that given by Fujii, but given the small- 
scale of the figures in [3] it was impossible to displa}^ the difference. Finally, it should be noted 
that Fujii estimated his numerical computations to be accurate to about one percent. Thus, 
the present solution is about as accurate as that of Fujii for Prandtl numbers of order unity. 
The Prandtl number restriction is necessary because the momentum and thermal plumes 
must be of comparable thickness for the above scheme to work. When Pr^>l, the momentum 
plume is much thicker than the thermal plume; while the reverse is true when Pr«Cl. 

The pressure variation within the plume may now be determined from the last of eqs 
(10). Using eqs (11), (13) and (14), the following expression is obtained: 



P= 



9o 



2j8(X-l) 



(1+%) 



-(X-l). 



1 / ay V 



(18) 



The function P(n) is plotted in figure 2 for Pr = 0.7 and 1.0. The profiles shown approximate 
the radial dependence of the pressure near the axis of symmetry far from the source. The 



I .70 



• 30 



0.90- 



P(rj) 



0.50 



0.10- 



-0.30 




0.00 



5.00 



10.00 15.00 



20.00 



25.00 



Figure 2. The function P^) describing the pressure profile inside the plume region. 
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initial drop in pressure as 77 decreases from infinit}^ is due to the inertial acceleration of the 
fluid as it enters the plume. The motion here is effectively irrotational, as will be shown in the 
following section. The pressure rise near the center is due to the buo}^ancy contribution, the 
first term in eq (18). A pressure rise near the axis of symmetry is also necessary on geometrical 
grounds as a means of turning the fluid from an inclined to a nearly vertical direction. A more 
complete discussion of the flow pattern and pressure distribution must await the calculation 
of these quantities in the region outside the plume where the inequality £^>?7 does not hold. 

4. The Uniformly Valid Pressure and Velocity Fields 

The determination of the pressure and velocity fields outside the plume can be separated 
from the calculation of the enthalpy. This separation is possible because the motion in this 
region is irrotrational. This may be verified by noting that the boundary layer solutions take 
the following form as 77— >co . 

G-goQrfJ^GWivr™^ (19) 

The last of eqs (19) follows from the fact that X>2 for the Pr range of interest. Substituting 
eqs (19) into eqs (8) and (9) the stream function and vorticity become: 

*=v (^J 2 £ 2 W (20) 

w=0. 

Thus, the asymptotic form of the solution valid in the plume is an irrotational flow. Now any 
irrotational flow is an exact solution of the Navier Stokes equations (1) provided the buoyancy 
terms may be ignored. Under these circumstances, the pressure is determined from Bernoulli's 
equation 

p-p m +lpv 2 =0. (21) 

In equation (21), the fact that u approaches zero outside the plume has been used. The 
velocities implied by eq (20) may be determined from eqs (9). Upon substitution into eq (21), 
the pressure becomes 



p-p-- 



oQ irafOT L__. (22) 

*- l - J (V+ 4 V) 



Comparing eq (22) with the last of eqs (8), and making use of eq (19), it is easily seen that 
eq (22) is consistent with the asymptotic plume solution when £^77. It is also consistent 
with the fact, noted in the previous section, that the inertial effects dominate buoyancy con- 
tributions to the function P(rj) for large rj. Finally, the velocities and pressures implied by 
eqs (20) and (22) satisfy the boundary conditions at infinity given in eq (2). Thus, the velocities 
given by eqs (9), (13), and (17) constitute a uniformly valid approximation to the velocity 
for £ 2 +l/4r? 2 ^>l (i.e., as \r\ approaches infinity in all directions). The corresponding uniform 
approximation for the pressure is readily seen to be 

gQ P(v) 
P—Poo=^r L ~ 



2 +k) 
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The stream function is of course determined by eq (8). 

The flow pattern is illustrated in figures 3 and 4 for Pr=0.7 and 1.0, respectively. These 
values are chosen because Pr=0.7 is appropriate for most common gases, while Pr=\ is an 
often used approximation in boundary layer calculations. The heat source is at the origin and 
the gravity vector points to the left. The bouyant plume occupies the region where the stream- 
lines are approximately parallel to the axis of S3^mmetry. The streamlines well away from the 
plume are parabolas of constant £. The transition region is observed to be roughly a parabola of 
constant 77. The streamlines are all asymptotically parallel, sufficiently far above the heat 
source. 
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Figure 3. The streamline pattern for several values of Figure 4. The streamline pattern for several values 
the dimensio riles s stream function (gQJph m v i ) l/2 'l'/v of the dimensionless stream function (gQ/ ph oo^) 1/2 ^/v 

and a Prandtl number of 0.7. and a Prandtl number of 1.0. 

The heat source is at the origin and the gravity vector points to the The heat source is at the origin and the gravity vector points to the 

left. left. 



The corresponding pressure distributions are shown in figures 5 and 6, respectively. The 
contours corresponding to p-p m ^>0 lie to the right (i.e. above) the heat source. These contours 
are all closed. They correspond to the positive portion of the P(fj) curve in figure 2. This high 
pressure region is dominated by buoyancy effects. The p=p a> curve is of course a parabola n = 
constant, the value of 77 being determined by the zero of P(r?). The remainder of the flow field 
is characterized by negative values of p-p m . This region is inertially dominated and as mentioned 
earlier, is nearly irrotational. The very rapid variation in pressure near the origin is probabty 
spurious, since the assumption that £ 2 +l/4r? 2 ^>l is certainly violated there. 

The overall flow pattern that emerges is thus determined by the relative magnitudes of the 
buoyancy and pressure forces. Inside the plume, the bouyancy dominates in the vertical mo- 
mentum balance, providing the force necessary to overcome the viscous drag and produce the 
upward motion. The lateral momentum balance is more complex, since the pressure acts first 
to draw the fluid into the plume, and then to retard and turn it toward the vertical. Outside 
the plume, the buo}^ancy rapidly becomes ineffective, since the lateral diffusion of heat is not 
sufficient to overcome its upward convection to the extent necessary to compete with the pres- 
sure forces. The pressure then acts to draw enough fluid into the plume to match the mass 
flux required at each station to carry the constant energy flux released by the source. 
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Figure 5. Contours of constant dimensionless over- Figure 6. Contours of constant dimensionless over- 
pressure hco(p — Poo)/gQ for a Prandtl number of pressure h. m (p — Pa>)/gQ for a Prandtl number of 
0.7. 1.0. 

The heat source is at the origin and the gravity vector points to the left. The heat source is at the origin and the gravity vector points to the left. 

The enthalpy outside the plume is not high enough to affect the motion to any significant 
extent. However, its determination is necessary to verify that the enthalpy distribution in the 
plume can be joined smoothly to a solution that vanishes as £ 2 +l/4r7 2 -^oo for all %jn. 
This solution cannot be obtained by the methods used to obtain ^ and p-p m . It is therefore 
necessary to discuss the sense in which these solutions satisfy the equations of motion, 
so that comparable approximations can be made for the enthalp}^. To this end, it should be 
noted that the plume region corresponds to the region $— > 0, where 4> is the polar angle measured 
from the positive z axis (see fig. 1). This may be seen by noting the transformations 

* 2 +i " 2 =(i^ 2 + 22 > 1/2 



2£_ 



= cot 



(!)■ 



Thus, the irrotational flow occupies the region > for distances (r 2 +z 2 ) 1/2 sufficiently far 
from the origin. The actual distance (r 2 +2 2 ) 1/2 required for £ + 1/4t7 2 ^>1 is dependent upon 
the parameter (gQlph ao v z ) li2 . The irrotational velocity and pressure were required to match 

the asymptotic forms of the plume solutions as $ and hence — vanish. They also were required 

to vanish at infinity. The same conditions will now be applied to the determination of the 
enthalpy. 



5. The Uniformly Valid Temperature Distribution 

The solution is based on the fact that outside the plume the stream function ^ may be 
approximated by eq (20). The energy equation, eq (7), then becomes: 



2^(co) _+- _ (^ _ j+ 2 _ (^ _ )= 
pv\pv 6 h m J 



(23) 
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H must vanish as £ 2 +l/4r7 2 -^oo. As <j>— > 0, the solution must match that given by the second of 
eqs (19), the asymptotic form of the solution in the plume. These boundary conditions may 
be written explicitly as: 

H=^( v y™^; |->0. (24) 

Once the solution satisfying eqs (23) and (24) is obtained, it will be combined with the solution 
obtained in section 3 in a straight-forward way to produce the desired result valid for all 
£ 2 +W»l. 

The solution to eq (23) may be obtained by Hankel Transform techniques. The Hankel 
Transform H(s,v) of H is defined by [13] 

/»00 

H(s, v )=j o dHJ (&)H<£, v) (25) 

r% oo 

//(f,i7) = J o dssJ o miT(s, v ). 

Multiplying eq (23) by the Bessel Function J (£s) and integrating over £ from zero to infinity 
yields: 

The solution of eq (26) satisfying the boundary condition at infinity is 

Tl=(r,s)-"K,(^jh( S y (27) 

Here K x is the Modified Bessel Function of order X[14]; and A is denned in eq (14). Thus, from 
eq (25) : 

II(Z,v) = ^dss J (m(s) (ns) Jx K, (|) (28) 

The function h(s) must be determined from the matching boundary condition given in the 
second of eqs (24). This boundary condition can be applied to eq (28) in a straight-forward 
manner if boundary conditions are known as rj — > for all £. Thus, letting r? — * in eq (28) 
with £ fixed by using the limiting form of K\(rjs/2) for small values of its argument: 

H(S,r,)=Qy\(\) (v)-*f°ds Joif) («) 1_a A(«) (29) 

Here, r denotes the Gamma Function [14]. The second of eqs (24) is now temporarily replaced 
by the following boundary condition, which is assumed to hold for all £ as 77 — > 0: 

ff =fq^ W~ 2X - (3°) 

Equating (29) and (30) : 

(3D 
A*(s)= s -2x h(s) 
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Equation (31) can be readily inverted making use of a special case of the Modified Weber- 
Schafheitlein integral [15]. 



/; 



JKo(os)«/o(^)s<fe=72X^ 2' ( 32 ) 

Thus 

h* = (2y-*Y^ K (s)- (33) 



Finally, using the inversion integral (25) : 

H(«, V)=^y »(») J " *««.(•) ('-J ^x (f ) Jo(«©. (34) 

Equation (34) constitutes an acceptable solution to the problem, since the matching 
condition, eq (24), is contained in the boundary condition actually satisfied. They become 
identical when £»1. However, there is a certain arbitrariness in eq (34). There are many 
functions which are regular for all £ and are proportional to (£)~ 2 for £M. An}^ such expression, 
when used in place of the factor (£ 2 +l) _1 in eq (30), would give a different solution for h*(s). 
It would be expected, however, that the behavior of h*(s) as s— >0 would be identical for all 
such choices, since the £~ 2 behavior as £— > °° must be preserved if the matching condition is to 
be satisfied. The arbitrariness can be removed by considering the limiting form of eq (34) when 
£ 2 +1/4t7 2 is large, but %/rj is arbitrary. This is accomplished by replacing the function K (s) by 
its limiting form for small s. The result is 



m > *?)= roT* G(co) L" dss log * (0* Kx (?) MH) ' 



(35) 



Now eq (35) is necessarily an exact solution of eq (23) since only h(&) has changed. It also 
satisfies the boundary conditions, eqs (24), as will be demonstrated below. 
Equation (35) can be rewritten in the form: 

J=f"ds(syK,(?£)J ( S !t) (36) 

Using the general form of the Modified Weber-Schafheitlein integral [15]: 

a=| (o+X+1) (37) 

The quantity 2 F'i is the Hypergeometric function [14]. Carrying out the differentiation, eq (36) 
becomes : 

ff-=^^[(iog(i)+ii*(i)+*(i+x)])(,.+«r»« 

+(.>-— Aj^M.-f))^., J 08) 

The quantity <f?(z) in eq (38) is the Digamma Function [14]; and the following result has been 
used: 

2F l (a, b, b;z)=(l-z)-. 
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In order to carry out the differentiation of 2 F\(a, b, 1, — 4£ 2 /»? 2 ), it is necessary to have 
an explicit representation of this function. The Barnes integral [14] is convenient, but can 
only be used to describe JFi{a, b, c, z) for c>6>0. Accordingly, the Gauss relations are used 
in the form: 

oJM, b, 1, z)=jZT z {[a-l-il-fyzhFiia, b, 2, z)-\-(2-a) 2 F l {a-l, b, 2, 2)}. (39) 

Thus, since b is in the neighborhood of 1 in the domain of interest, the Barnes integral can be 
used on the right-hand side of eq (39) to obtain: 

2 F,(a, b, 1, z)=(i- 2 )- 1 {r(6)r(2-6)}- 1 . 

j[(o-l)-(l-6)«] £ 1<fe (rZi) >_l (l-^)- a +(2-«) £ dx {jZ^)*~ 1 (l-*2)- ( - 1> }- (40) 
Thus: 
§j*Fi(fl,b, l,2)| p . I+x = [2(l-2)]- 1 

x{(X 2 )- 1 (l + 2)[(l-3)- X -l]-[(X-l)2]- 1 [(l-s)- (X - 1) -l] 

-(X2)-'{[(l-2)- X -l] + X(l-2)- X log (1-2)} 

+ [(X-l)2]- I {[(l-3)- <X - 1) -l] + (X-l)(l-2)- <X - 1) log(l-2)} 

+ \Cdx log (V^) (l-^2)- a+1, + (l-X) f efe log (jz^) (l-2^)- x i • (41) 
Inserting eq (41) into eq (38) with 2= — 4£ 2 /?? 2 ; the following expression is obtained: 

H-f=^{2Xlog(j)(l+«-H/(0 

/(t) = -A[*(l)+¥(X+l)] (l + < 2 )^ + ^!zl) [l_(l + ««)-*] 

+x[(x-i)< 2 ]- 1 [i-(i+« 2 )- (X - 1> ] 

+ | {l-(l + i 2 )- x -X(l + ( 2 )- x log (1 + i 2 )} 

+2X[(X-1)« 2 ]- 1 {(1 + < 2 )" <X_1) -1 + (X-1) (l + ^ 2 )- w - 1 » log (l + t 2 )} 

-?|a+«-^|[(i+^-i]-X ,, |[i-(i+*)- k ]} 

+ Aj(l+^)-(x-i)J^[(l +a; )x- I _i]_JJ J ^ [ i_ ( i +a; )-(x-i) 1 j 

t =f=cot(|) (42) 

The demonstration that the boundary conditions, eqs (24) are satisfied can be made using 
eq (42). Since t depends only on the polar angle <j>, H obviously vanishes as £ 2 4-l/4 j? 2 ^co for 
any non-zero value of $ (finite t). Furthermore, f(t) has the limiting behavior 

/«) = -A[*(l)+*(A+l)] + 6>(( 2 ); t^O 
j(t) = l + 0(t- 2 ); «->•. 
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As the plume region is approached (<£— >0; t— ■> °° ), the matching boundary condition is recovered. 



Finally, a uniform approximation to H(%, rj) valid for all <f> may be constructed. Note that 
the quantity 6?(°°) (^)~ 2X appearing in eq (42) is the asymptotic form of the plume enthalpy 
function G(rj). Now consider the form 



H(€,n)= 7 fl ' ( ? ) x fx log lYjY+ll (l+* 2 )- x +/«)] 



(43) 



Applying the limit £-><» with 77 fixed to eq (43), the form for h— h m given by eq (8) is recovered. 
The limit £+-- r/ 2 -^ 00 with <p fixed yields eq (42). The replacement of 2log (77/4) by log [(r?/4) 2 +l] 
is motivated by the desire to have the uniform expansion be an analytic function of n 2 . 

The function /(£) is plotted as a function of polar angle <j> (eq (42)) in figure 7. The function 
is nearly constant over the region occupied by the plume, but varies rapidly elsewhere. Figures 
8 and 9 illustrate the enthalpy distribution, showing several level curves of the function H(£, rj) 
given by eq (43). Again, the detailed structure near the origin cannot be taken seriously, due 
to the smallness of £ 2 + l/4r? 2 . The figure does illustrate, however, the rapidity with which 
the excess enthalpy decreases away from the plume region near the positive z axis. 



/>r = IO 



f(<£) 




</> 



Figure 7. The function {(<f>) appearing in eq (43) plotted as a function of polar angle measured in degrees. 

<f>=0 corresponds to the positive axis (the plume centerline). 
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Figure 8. Contours of constant dimensionless enthalpy Figure 9. Contours of constant dimensionless enthalpy 

rise pKpj^hco/gQ) 1 / 2 ^ — tO/Q for a Prandtl number rise pv(pifih m /gQ l / 2 (h — h m )/Q for a Prandtl number 

of 0.7. of 1.0. 

The heat source is at the origin and the gravity vector points to the left. The heat source is at the origin and the gravity vector points to the left 



6. Discussion 

The present results shed some light on the use of the Oseen approximation in buoyancy 
induced flow problems. In the present context, the "Oseen approximation" is taken to mean 
the replacement of the convection velocity b}^ a substitute quantity whose choice simplifies 
the calculation of the flow quantities of interest [16]. As noted in the introduction, Fendell 
[11] and Wesseling [12] use this technique to study convection induced b}^ a localized heat 
source. Wesseling assumes the heat source to be at rest in a uniform stream produced by some 
external agent, whereas FendelFs motion is induced entirely by the heat source. 

First consider FendelPs analysis. Here, the convection velocity is assumed to be the center- 
line plume velocity, with the result that the velocity and temperature perturbations decay 
exponentially with distance away from .the plume. The actual convection velocity outside the 
plume is given by the second of equations (9) with F(rj)=F(co). This velocity itself decays 
slowly in all directions, and is much less effective at inhibiting the diffusion of heat away 
from the plume than a constant velocity parallel to the plume would be. 

The analysis contained in the previous section ma}^ be thought of as an Oseen linearization 
with the irrotational flow as the convection velocity. The irrotational flow is dominant in this 
region, as discussed in section 4. However, there remains a residual vorticity containing motion 
outside the plume, analogous to the enthalpy distribution calculated in section 5. This motion 
could in principal be calculated b}^ applying an Oseen linearization to the vorticity equation 
with the convection velocity described above. In the author's view, such a calculation would 
miss the main point; that the dominant motion outside the plume is irrotational. This irrota- 
tional motion is not known a priori and must be determined from a knowledge of the plume 
solution, including the pressure distribution. Once it is known, the homogeneous solution to 
the vorticity equation, corresponding to co = 0, is known. The solution co = thus implies the 
non-trivial flow discussed in section 4. It is the existence of the irrotational flow which dis- 
tinguishes the solution o> = of the vorticity equation from the trivial solution h—h m =0 of 
the energy equation. This is why the small quantity h—h^ must be calculated outside the 
plume, while the corresponding vorticity distribution is unnecessary unless desired for its 
own sake. 
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The preexisting uniform flow in the absence of any heat source is used as the convection 
velocity in Wesseling's analysis. The consistency of his result then requires the buoyancy in- 
duced velocities to be small compared with the uniform flow everywhere away from the source. 
Two cases are considered; a point source and a line source. Wesseling finds that while the point 
source has an acceptable solution, the line source does not. He remarks "it remains to be seen 
whether the Stokes-type paradox exhibited above is caused by the Oseen linearization or 
whether it is a property of the full Boussinesq equations." Although it is not possible to settle 
this question here, the following facts are suggestive. The present analysis together with pre- 
vious work demonstrate that the thermally induced velocity in a medium at rest is bounded 
everywhere away from a point source. On the other hand, the velocity of a plume generated by 
a line source in a medium at rest increases with the % power of the distance above the source. 
This must induce a velocity outside the plume which is ultimately larger than any uniform 
velocity. Thus, the thermally induced velocity probably dominates the motion outside the 
plume in a moving medium in two dimensions. This velocity cannot be assumed a priori as was 
done by Wesseling, but must be determined. If detailed analysis bears out the above hypothesis, 
then it is the Oseen linearization which must be reconsidered in two dimensions rather than the 
Boussinesq equations. 

The author is indebted to Margaret Harkleroad for the numerical computations and the 
figures appearing in this work. 
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